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Abstract: 

In the final stages of stellar evolution low- to intermediate-mass stars lose their envelope in increasingly 
massive stellar winds. Such winds affect the interstellar medium and the galactic chemical evolution 
as well as the circumstellar envelope where planetary nebulae form subsequently. Characteristics of 
this mass loss depend on both stellar properties and properties of gas and dust in the wind formation 
region. In this paper we present an approach towards studies of mass loss using both observations and 
models, focusing on the stage where the stellar envelope is nearly empty of mass. In a recent study we 
measure the mass-loss evolution, and other properties, of four planetary nebulae in the Galactic Disk. 
Specifically we use the method of integral field spectroscopy on faint halos, which are found outside 
the much brighter central parts of a planetary nebula. We begin with a brief comparison between our 
and other observational methods to determine mass-loss rates in order to illustrate how they differ and 
complement each other. An advantage of our method is that it measures the gas component directly 
requiring no assumptions of properties of dust in the wind. Thereafter we present our observational 
approach in more detail in terms of its validity and its assumptions. In the second part of this paper 
we discuss capabilities and assumptions of current models of stellar winds. We propose and discuss 
improvements to such models that will allow meaningful comparisons with our observations. Currently 
the physically most complete models include too little mass in the model domain to permit a formation 
of winds with as high mass-loss rates as our observations show. 
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1 Introduction 

Low- to intermediate-mass stars, with a main sequence 
mass of 0.8 Mq <, M <, 8Mq, lose most of their en- 
velopes as they ascend the thermally pulsing asymp- 
totic giant branch (TP- AGB). The increasing rate of 
mass loss dominates time scales of stellar evolution as 
the stellar wind changes the star more rapidly than nu- 
clear burning in the stellar core does. The matter of 
the stellar wind is enriched with nuclearly synthesized 
elements, which are dredged up from the stellar inte- 
rior during the evolution. These winds are also sites of 
strong dust formation and determine properties of the 
physical environment where planetary nebulae (PNe) 
form subsequently. Because most stars pass through 
this phase mass loss from AGB stars contribute sig- 
nificantly to the chemical evolution of the interstel- 
lar medium and the galaxy. Recent reviews of mass 
loss from cool stars are given by, e. g.. IWillson 
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iHofneil l|2005r i. IZiilstral (|2006l '). and Ivan Loo: 
Also see Zijlstra (these proceedings). The review of 
iHSfncr (2009) focuses on the role of dust. 

Stellar winds of AGB stars form in a complex en- 
vironment of many variables, here is a brief and much 
simplified description of the formation process. These 
stars tend to pulsate as long period variables, with 
time scales of several hundred days. The pulsations 
form in the outer part of a convective envelope, that 
surrounds a bright giant star where hydrogen and he- 



of about lO'' 



-10* yr- 



The stellar atmosphere is cool 



and allows formation of large amounts of sub-micron- 
(to micron-) sized dust, which is either carbon-rich 
(C-stars) or based on other elements (M-stars and S- 
stars). The dust forms under non-equilibrium condi- 
tions at similar time scales as the pulsations, and dust 
in many circumstances provides the main driving force 
for the wind due to its large opacity; the dust-driven 
wind forms when individual dust particles drag the 
gas along by collisions. The wind shapes an extended 
stellar atmosphere, and further out a circumstellar en- 
velope where the radial (r) density gradient would be 
proportional to at a constant mass-loss history. 

When the central star is not concealed by a mas- 
sive envelope anymore it heats up and can, under the 
right circumstances, ionize the matter of the AGB 
wind. The shape of the then formed ionized (plan- 
etary) nebula shows a variety of geome tries (see e.g. 
iBalick fc Fraii^ l2002l : ISahai et al.l l2007D . In view of 
PNe with a spherical or elliptical geometry they con- 
sist of a central part with a relatively bright rim and a 
fainter attached shell; both are built-up in an interplay 
between ph oto-ionization and wind interaction (see, 
for example. [ Fran k et al.||l990|: [Marten fc Schonbernerl 
ll99ll : lFrankll994l : lMellemall994l : lPerinotto et al .1120041 ). 
The matter in the central parts of these objects has 
been restructured during their formation and this mass 
can therefore not be used to measure properties of the 
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last period of more spherically symmetric mass loss on 
the AGB. As the PN evolves the region outside the PN 
proper is eventually ionized, and a halo forms (also see 
Sect. [2]) • Such halos are relatively unaffected by pro- 
cesses of the formation of the PN, and therefore still 
contain directly accessible inf ormation of the previous 
AGB mass loss. ISokerl 1120061 ) discusses the transition 
from more spherically symmetric outflow on the AGB 
to axi-symmetric outflow that is often observed in later 
stages. 

Measurements of mass loss from AGB stars are 
made by four different methods, using spectroscopy, 
imaging, or both. Observations, which cover and spa- 
tially resolve the circumstellar envelope, readily pro- 
vide information on the mass-loss evolution through- 
out the envelope. 

In one spectroscopic method radio emission lines 
of carbon monoxide are compared with models of ra- 
diative transfer in the circumstellar envelope. Mass- 
loss rates ar e here by calculated using best-estimate 
fits; ISchoierl (|2007[ ) reviews this method. The relia- 
bility of the outcome of this method, that measures 
the gas-com p onent in the w ind, is discussed by e.g. 
ISahai l|l99d ). iKastnerl ||1992D , and lSchoier fc OlofssonI 
l| 200 il l . This approach works best with mass- loss rates 
M^lO~^M0yr~^ (at higher values the relevant lines 
become saturated: iRamstedt et al]l2008l l. 

In another spectroscopic method infrared measure- 
ments of dust are used to match the observed spec- 
tral energy distribution with mode ls of distr i butio n 
of dust; this method is reviewed bv Ivan Loon! (|2007ri . 
Different models are used to prescripe input mass-loss 
rates, rangin g from constant valu es to mass-loss histo- 
ries (see e.g. lGroenewegeiilll995ll . In order to fit the 
spectral energy distribution of some stars in the Mag- 
ellanic Clouds it seems necessary to vary the model 
input mass-loss rat e by a factor of a few ove r time pe- 
riods of 10^-10^ yr (jGroenewegen et al.ll2007l . although 
as the authors comment this should be a less likely sce- 
nario). Both spectroscopic methods require additional 
assumptions of properties of the mass-loss rate, viz. 
the expansion velocity, and the dust-to-gas mass ratio. 

iMarengol (|2009l ') reviews imaging methods for mea- 
suring the mass loss, the chemistry, etc., using filters 
covering ultra-violet to infrared wavelengths. 

In a fourth method lSandin et al.l (|2008l ) study four 
halos of Galactic Disk PNe and measure the mass-loss 
evolution and other properties after the last thermal 
pulse near the tip of the AGB. Our new approach 
uses the novel method of integral field spectroscopy 
(IFS), that in principle is a low spatial re solution imag- 
ing rn ethod. Based on the approach of iPlait fc Sokerl 
(| 19901 ) this method allows the radial gas density struc- 
ture - i.e. the mass-loss evolution - to be measured 
directly. No additional fitting procedure is necessary 
using supplementary models, as in the other two spec- 
troscopic methods. Lacking better data, however, a 
constant value is currently used with the expansion 
velocity. This method, furthermore, only works when 
observing the final stages of AGB mass loss since the 
matter only becomes ionized during the PN phase. It 
obviously also only works with stars which form a PN 
(and later also a halo) . The authors find that the mass- 



loss evolution shows increased values by a factor of 
about 4-7 during 10* yr for all four objects. 

iRenzinil (Il98ll . also see lPeimbertlll98]l ) introduces 
the term superwind in order to emphasize that stars 
ascending the AGB reach higher mass-loss rates than 
stars on the precedent red giant branch. Using stellar- 
evolution calculations including an empirical mass-loss 
relation (that depends on the puls ation period, but not 
the st ellar mass or the luminosity) IVassiliadis fc WoodI 
1.1993!) find that mass-loss rates reach the hig hest val- 
ues during the final thermal pulses on the AGB. lBlockeri 
( 1995 , who also provides a summary of previously de- 
rived mass-loss relations) uses a mass-loss relation for 
the AGB, that is based on the red-giant branch relation 
of Rcimcrs ( 1975, alth ough account ing for effects of the 
dynamical models of iBowenI Il988ll , and thereby also 
accounts for increasing mass-loss rates i n his stellar- 
evolution calculations. ISchrode~ et all (|1999D reach 
the same conclusion, using a mass-loss relation that 
is based on models of dust-dr iv en winds (adopt i ng the 
relation of lArndt et~allll997r ). Ivan Loon et al.l (| 19991 . 
120051 ). furthermore, come to a similar conclusion with 
their observationa l st udies of dust-enshroude d stars. 
iHeske et all (|1990| ) and lJusttanont et al.l (|1996D present 
additional results which support the above findings. 
Using timo-dcpendondont photoionization models of 
PNel^chonberncr ct al. (2005), moreover, find that the 
AGB wind must increase in strength during the final 
evolution. Models assuming a stationary AGB wind 
fail to explain virtually all observed PNe. 

Due to the small number of assumptions our obser- 
vational method provides a powerful approach to study 
the decisive stage of stellar evolution with strong mass 
loss at the tip of the AGB. In this paper we present the 
key concepts of our method, including suggestions on 
how to use it with additional objects. We also discuss 
required improvements to models of stellar winds to 
be used when interpreting the outcome of our observa- 
tions. We first describe properties and observables of 
the objects of our observational method, AGB halos of 
PNe, in Section [2| Thereafter we discuss assumptions 
and applicability of current models of stellar winds in 
Section [3| 

2 About AGB halos of PNe 

The definition of what a PN halo is comes from photo- 
ionization models of PNe, which have a spherical or 
elliptical shape. The halo is identifi ed as the ionized re- 
gion outside the central PN (see e.g. lSteffen fc Schonbernerl 
120031 : IViUaver et alll2002al lbl) . It appears a few thou- 
sand years after the central PN becomes visible. The 
halo is also less affected by PN shaping processes than 
the central region is, and therefore allows direct mea- 
surements of properties of the AGB wind. However, 
in general the least affected region of mass loss in the 
halo is the inner part. In the outer halo it is expected 
that wind interaction, between periods of (the last) 
thermal pulse(s) and intermediate periods of mass loss, 
has modified the fiow. The interaction between the PN 
and the inters tellar medium (ISM) is studied by e.g. 
IVUlaver et all (|200a ) and .Wareing et al.. (,2007i ). They 
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find that the nebula can become asymmetric when the 
sphericaUy symmetric structure of an AGB wind moves 
through the ISM. For s uch objects (which fr action is 
about 20% of all objects. IWareing et al.l[2007l l it is im- 
portant to measure the mass loss in the wake of the 
interaction. 

2.1 How to differ between AGB ha- 
los and other structures 

The term halo has in the past oft en been used to indi- 
cate different types of structures. ICorradi et al.l (|2003l . 
cf. Section 3.2 therein) require four criteria to be met 
in order to positively identify the presence of an AGB 
halo. We use the same criteria, but rephrase them 
slightly as follows: 

• The halo is limb-bright ened, or de t ached accord- 
ing to the definition of IChu et all iIIQSTT ). 

• The surface brightness (SB) of the halo must 
be of the order lO'^-lO* weaker than the peak 
emission on the rim in the inner nebula, or even 
weaker. The intensity of the halo decreases with 
increasing radius, and halos of some objects show 
this full intensity range. 

• The radius of the halo is at least twice as large 
as the inner rim radius. 

• Objects which appear to be evolved and have a 
low-luminosity central star together with an av- 
erage halo-SB of about 10~^ of the peak inten- 
sity are candidates to have a recombin ation halo 
fcf. lTvlenda|[l986l : ICorradi et al.ll2000l ). Such re- 
gions are not fully ionized anymore, and in com- 
parison to younger PNe the matter in the halo 
has been restructured during the evolution of 
the PN. These objects provide little information 
about the AGB wind. 

The bright central parts of a PN may infiuence mea- 
surements of t he halo around it through i nstrumental 
scatte red lig ht. iM iddlema ss et al.l (|l989al ). ICorradi et al.l 
(|2003l '). and ISandin et al.l (|20oi ) present methods to 
dis tinguish halos f rom in strumental scattered light. 

ICorradi et aD (|2003h find 25 Galactic Disk AGB 
halos, which are either circular, slightly elliptical, or 
highly asymmetric. This is about 60% of the sam- 
le of elliptical objects in their studv. [ Sandin et al.l 
20081 ') find two additional halos; ICorradi et al.l (|2003l ') 
mark one of th e se ob jects, IC 3568 to not have a halo. 
iReid fc Parkeil l|2006l . SMP27) and Cohen (these pro- 
ceedings) identify two additional AGB halos in the 
Magellanic clouds. Although relatively few AGB ha- 
los are known it appears that they are a common phe- 
nomenon, where their detection is mainly limited by 
their brightness. 

2.2 About rings found around AGB 
stars, post-AGB stars and PNe 

In many cases images of AGB stars, post-AGB stars 
and PNe reveal multiple ring-shaped shells (also called 
arcs) around the main object. These shells are sep- 
arated by time periods of about 10^-10^ yr, see e.g. 



ISahai et al.l (ll998l.ll999fl.lKwok et al.l (119981) . ISu et al l 

(1998). Mauron fc Huggins (1999. '200tf. 2006"). Bond! 
mOO). Hrivnak et al. (2001), Terzian_fc^Hajiaa ( 200cj ), 
and lCorradi et al.l (|2003. 2004). The origin of the rings 
is not understood, but numerous processes have been 
suggested to explain their fo rmation. For example, 
iMastrodemos fc MorrisI (Il999|) argue that the rings re- 
sult from binary interaction. Isimis et al.l ([2001) ad vo- 



catc an origin from a drift-induced instabilitv. lZiilstra fc Bedding 
(2002 | ) suggest a pulsation al instability as the origin, 
and Ivan Horn et al.l (|2003l ') argue for a nuclear- burning 
instability. A more complete discussion of these and 
additional suggested mechanisms is given bv lSokeil l|2002l . 
[2006^). Our method to measure mass-loss rates using 

IFS results in data with too low spatial resolution, in 

gener al, to find these rings (see, however. iMonreal-Ibero et al.l 
l2005l . for NGC3242). Although the rings arise due to 
stellar processes, which are important to understand, 
their infiuence on measurements of the mass-loss evo- 
lution currently appears minute. In order to measure 
and evaluate density differences of matter, which is 
present in the rings and between the rings, it is nec- 
essary to further increase the spatial resolution of the 
mass-loss measurements. 

2.3 Observables in the halo for the 
optical wavelength range 

The low surface brightness of AGB halos implies a 
lower number of observable emission lines than in the 
bright central PN. Visible emission lines in the opti- 
cal wavelength range, which span most of the radial 
extent of the halo, include the Balmer lines Ha and 
H/3, the two bright oxygen lines [O ill] AA4959, 5007, 
and often also neon [Ne ill] A3869, 3968 ([Neiii] A3968 
may be blended by He). In the innermost part of 
the halo the number of visible lines is larger. De- 
pending on the density and the ionization structure 
some of the density-sensitive lines [O ll] AA3726, 3729, 
[Ariv] AA4711, 4740, and [Sii]AA6717, 6731, are nor- 
mally visible. The weak temperature-sensitive aureole 
line [O III] A4363 is mostly visible, although depend- 
ing on the site it may be blended by a much stronger 
telluric mercury line Hg A4358. Additional emission 
lines, that are often visible, include H7, H5, Hel A4471, 
HeiiA4686, and possibly [Nil] AA6548, 6583. 

The low density throughout the halo requires it to 
be measured using recombination lines. Hereby the 
strongest Balmer lines can be used to calculate the 
electron density in the halo (assuming a distance to the 
object) and thereafter the gas density and the mass- 
loss rate. The line ratios of [0 11], [Ariv], and [S 11] can 
be used to scale the density structure using values mea- 
sured near the inner boundary of the halo. In our study 
it turned out that the recombination-line density struc- 
ture matches the line-ratio densities only for one object 
assuming a proper distance (viz. NGC6826). For the 
remaining three objects it is, for some reason, neces- 
sary to use much shorter distances when calculating a 
recombination-line density structure that matches the 
line-ratio density measurements. It therefore appears, 
that because of this scaling it is, in general, not possi- 
ble to determine an object distance. 
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The oxygen lines are normally used to calculate 
the electron temperature in both the central PN and 
in the halo; the temperature-sensitive nitrogen line 
is too weak in the halo. Detailed radial structures 
of the temperature in several cases reveal hot halos 
(ISandin et al.ll2008|. such tendencies are also identified 
bylMiddlemass et al."l989b','l99l';'Manchado fc Pottasch| 
1989 : Monrcal-Iboro ct al. 2005; Krabbo & Copotti 20oj). 
MartenI (jl993 ) finds a simple explanation to the high 
temperatures with non-equilibrium models which ac- 
count for hydrodynamic effects. With such models the 
ionization front, that is responsible for the high tem- 
perature in the halo, moves through the envelope in a 
short time period (typically lO'^-lO'^ years long) . These 
models suggest that high temperatures in the halo are 
a fairly com mon, although trans ient phenomenon (cf. 
Sect. 6.1 in ISandin et al.l l2008l ). Accurate electron 
temperatures are necessary in order to measure abun- 
dances in the halo correctly - although, unless addi- 
tional emission lines at either shorter or longer wave- 
length ranges are used it is difficult to get such values 
for many elements. 

Th e four mass-loss evolution structures that lSandin et 
(|2008l ') measure assume a constant outflow velocity of 
lOkms"^. Mass-loss estimates will be more accurate if 
a radial outflow velocity structure could be measured 
also in the halo. The resulting line broadening, due 
to the expected velocities, is modera te. Nevertheless 
the approach of lGuerrero et al] (Il998l ) could be useful 
to make such measurements. Note that radial density 
and temperature structures, that are both provided by 
our method, are required in order to deconvolve kine- 
matic measurements when determining a radial veloc- 
ity structure. 

3 Modeling AGB mass loss 

A detailed theory of formation of stellar winds from 
cool stars should be able to reproduce observational 
findings. Because there are so many relevant physical 
processes in the AGB wind - where many processes are 
not well understood - a complete and self-consistent 
modeling of the AGB wind is (still) difficult to make. 

A complete model of mass loss on the AGB should 
be three-dimensional. It should cover a radial region 
extending from the core (and possibly even include the 
core) out into, and including, the wind formation re- 
gion, say, out to about 50 stellar radii. It should, more- 
over, properly handle the convective mantle between 
the core and the atmosphere, and thereby be able to 
describe the formation of stellar pulsations. It should 
handle dust formation in different chemistries as well 
as the gas-dust interaction. Radiative transfer should 
be done in a frequency-dependent approach account- 
ing for possible non-local thermal equilibrium. And 
the model should work for a significant period of the 
evolution of the AGB star. 

3.1 State of current wind models 

The formidable problem of stellar wind formation can 
only be solved by breaking it down into manageable 
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pieces, often with large - and necessary - simplifica- 
tions. Here follows a general list of such key pieces, 
which have contributed to the development of the phys- 
ical part of models of mass loss so far. The first time- 
dependent stellar wind models that ca n acco u nt for 
effe cts of pulsatio ns are presented by IWoodI (|l979l ') 
and iBowenI (| 19881 . who also, because of a lack of the 
necessary data at the time, int roduces a const a nt gas 
opa city, that stil l is in use). lFleische~ et al.l (|1992| ) 
and lHofner et all (Il995l ) present time-dependent wind 
models, whic h use the time - depen dent dust formation 
description of lGaueer et al.l (Il990l. for C-ric h dust; also 
see the references t herein). iHofner et al.l (1998) and 
iHelling et al.l (|2000l l introduce gray m olecular opaci- 
ties, w hich replace the constant opacity. iGail fc Sedlmavii 
(!l999|) extend the dust formation descripti on to work 
in M-type stars. ISandin fc Hofne^ (|2003l ) introduce 
the physics and role of gas-dust drift - although the 
physical outcome is somewhat affected by numerical is- 
sues (which are now corrected, see below) . iHof ner et al.l 
(2003) introduce freque ncy-dependent opa cities for both 
the gas and the dust. IJeong et al.| (|20 03|) study dust 
aljormation in M-type stars. [Sandin fc Hofnerl (|2004l ') 
find that gas-du st drift increases the amou nt of formed 
dust. Moreover, IWoitke fc Niccolinil (|2005l ) and lWoitkd 
(2006) study the outcome of two-dimensional wind mod- 
els. .Sandin (20081 improves the numerical accuracy 
of both drift and non-drift models (cf. Section 13. 2|) . 
iFrevtag fc Hofneil ([2008 1 present a first step tow ards 
three-dimensional wind models, an d iHofned (|2008l l flnds 
that a frequency-dependent approach and larger (micron- 
sized) grains may provide the necessary force to drive 
a wind off M-type stars. 

Self-consistent models of AGB mass loss should 
rely on a small number of input parameters - the 
mass-loss rate and the expansion velocity must not be 
such paramet ers. Current models, as described by e.g. 
'Sandin' (2008!), require the following four parameters to 
be specified: the stellar luminosity, the effective tem- 
perature, the composition (in this case the C/O-ratio), 
and the stellar mass. In one-dimensional models it is 
also necessary to describe the stellar pulsations using 
either a description for turbulent energy (i.e. convec- 
tion) or a moving boundary (a so-called piston). This 
adds about 2 parameters (cf. Section [3. 3|l . 

There are, however, more assumptions built into 
the models. The gas is described by the equation of 
state, assuming an ideal gas with a mean molecular 
weight. The interaction between the gas and the radia- 
tion field is described by opacities, which are calculated 
in advance for the necessary atoms and molecules. Al- 
though the gas opacity does not provide a strong enough 
force to drive a wind it is a decisive component that 
determines the density gradient in the atmosphere - 
and thereby influences the wind formation indirectly. 
In a sense, using a constant gas opacity all other de- 
tailed physics becomes of secondary importance. A 
number of simplified mass-loss rate formulae are in use, 
which are based on mass-loss models us ing a constant 
gas opacity; see e.g. lArndt et al. l (1 1997) , who se mass- 
loss relation is used by Schrod er et al.lll999l i n thei r 
stell ar-evolution ca lculations) , IWachter et al.l 12002") , 
and IWachter et al.l (2008 '). Because models using a 
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3.3 Modeling mass-loss formation 
near the tip of the AGB 



constant gas opacity achieve unreahstic density struc- 
tures, these formulae cannot be expected to provide 
quantitatively useful estimates of mass-loss rat es. Based 
on thi s argument alone the mass-loss relation of lMattsson et 
(|2009l '). that is provided in the form of an interpolat- 
ing computer program, seems more appropriate as it 
is inferred using frequency-dependent opacities. 

The dust component is, furthermore, likewise de- 
scribed by properties of the dust grains and through 
the du st opacity; in their parameter studv lAndersen et al.l 
(|2003l l find that wind models are highly sensitive to 
how such dust properties are used. The term describ- 
ing the interaction between gas and dust, the drag 
force, is rather insensitive to which type of colli sions 
are used (specular or diffusive) . ISimis et all (|200ll ) find 
that a gas-dust drift-instability is the cause of the rings 
observed around many (former) AGB stars (cf. Sec- 
tion [53}. Their models are, however, very difficult to 
reproduce as they are based on a descripti on using sev - 
eral simplifications (this is also noted bv lSokeill2002l l. 
The inner boundary is hydrostatic, thus pulsations are 
not considered. They treat radiative transfer in a sim- 
plified approach assuming an optically thin envelope, 
and use only one parameter setup of rather unrealistic 
properties. An extended period of mass loss is mod- 
eled by replenishing mass at the inner boundary (by 
keeping the gas density fixed). Instead of reproducing 
this setup we propose to make the required physical en- 
hancements to our models (see below) and thereafter 
see if and how a similar drift-instability appears. 



^order to model stellar winds near the tip of the AGB 
it is necessary to enhance current wind models. Nor- 
mally a radial domain is used that extends from above 
the region where pulsations form (at about 0.8-0.9 stel- 
lar radii) to 10-60 stellar radii. The amount of mass 
present in the model within this domain (about 10~^- 
1O~^M0, using gray or frequency-dependent opaci- 
ties, respectiveljij) is not enough to calculate mass-loss 
rates of the order ^10~^ M0yr~^ with any reliability. 
For example, it normally takes about 10^ years of evo- 
lution to start a dynamic model from a hydrostatic 
initial model. 

By including more of the central parts of the star 
in the model domain it becomes possible to model 
also the highest mass-loss rates - and thereby pro- 
vide a mean to understand our observations at the tip 
of the AGB. The inclusion of the convective stellar 

mantle in the model requires the use of eit her a three- 

dime nsional approach (such as presented bv lFrevtag fc Hofnerl 

200811 or a description for turbulent energy (see e.g. 

Kuhfudl98d:lGehmevr fc Winkl eJl99l:ICanuto fc Dubovikovl 
1998l : IWuchterl fc Feuchr ingcr 1 9981: iFeuchtingeJ 19991 ) . 
which implies a number of additional free parameters. 
It should, however, be possible to set many of these 
new parameters to default values. ISokeil (|2006) also 
suggests an approach where more of the entire star 
is modeled in order to fully understand the formation 
and role of the pulsations. Finally, it may also turn 
out that dust cannot itself drive the strongest winds 
observed, and in this case it may be necessary to find 
another, complementary, driving mechanism. 



3.2 Effects on models of improved 
numerical accuracy 

The numerical description in AGB star wind models 
is explored in detail by Sandin (2008). He introduces 
and studies several improvements to the numerical ap- 
proach. The results show that most one-dimensional 
wind structures are steady, with small periodic fluc- 
tuations. In comparison a majority of previous wind 
model structures are mildly chaotic. The formation of 
dusty shells in the wind, that were the topic of many 
papers up to now, do not form in the new models. Al- 
though, models which still use a constant gas opacity 
instead of molecular opacities continue to form chaotic 
structures; probably as a consequence of the fewer de- 
grees of freedom which are available in such models. 

Models which are calculated using gas-dust drift 
require a more careful numerical treatment than non- 
drift models do. It turns out that the dust is more 
tightly coupled to the radiation field than it is to the 
gas component. The dust thereby pulsates markedly 
along with the modulated radiation field. It should be 
remembered that in comparison to real (three-dimensi- 
onal) stars these models are one-dimensional. Effects, 
such as convection and formation of clumps, which 
could cause an inhomogeneous or even chaotic flow, 
are not considered in these models. 
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